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Abstract 

This review focuses on the processes that energize and trigger M- and X-class so- 
lar flares and associated flux-rope destabilizations. Numerical modeling of specific 
solar regions is hampered by uncertain coronal-field reconstructions and by poorly 
understood magnetic reconnection; these limitations result in uncertain estimates 
of field topology, energy, and helicity. The primary advances in understanding field 
destabilizations therefore come from the combination of generic numerical experi- 
ments with interpretation of sets of observations. These suggest a critical role for 
the emergence of twisted flux ropes into pre-existing strong field for many, if not all, 
of the active regions that produce M- or X-class flares. The flux and internal twist 
of the emerging ropes appear to play as important a role in determining whether an 
eruption will develop predominantly as flare, confined eruption, or CME, as do the 
properties of the embedding field. Based on reviewed literature, I outline a scenario 
for major flares and eruptions that combines flux-rope emergence, mass draining, 
near-surface reconnection, and the interaction with the surrounding field. Whether 
deterministic forecasting is in principle possible remains to be seen: to date no re- 
liable such forecasts can be made. Large-sample studies based on long-duration, 
comprehensive observations of active regions from their emergence through their 
flaring phase are needed to help us better understand these complex phenomena. 
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1 Introduction 



We have known of the phenomenon called 'solar flare' ever since the first 
observation by Carrington and Hodgson in 1859. We have learned a great 
deal in the century and a half that have passed since, but these basic ques- 
tions remain with us: What powers a solar eruption? What triggers it? These 
questions continue to occupy us as we bring ever-increasing instrumental and 
computational capabilities to bear on the problem. At the time of this writing, 
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in the summer of 2008, a search with the Astrophysics Data System yielded 
over 11,000 refereed publications with the words 'flare' and 'Sun' or 'solar' 
in the abstract. Even the limited focus of this paper on the energy source 
and triggering process of active-region flares and their associated eruptions is 
associated with thousands of studies. 

This review focuses on active-region M- and X-class flares and the flux-rope 
destabilizations associated with them. I present my summary first, followed by 
a synthesis scenario for the driving of field eruptions. These two sections are 
meant to give the reader a context to assess the supporting evidence that is 
discussed in detail, with references, in subsequent sections. Readers who would 
rather see the more traditional review in which the discussions and conclusions 
follow the arguments may skip Sect. 2 initially, returning to it just prior to 
Sect. 13. 

A selection of characteristic flares and eruptions discussed in the text is 
compiled in Table 1. These flares are all well-observed examples that illustrate 
the findings described in the indicated parts of the text, and many are the basis 
of studies that are discussed in this review. 



2 First, the conclusions: 

2.1 Assessing key findings and inferences 

One would like to have an understanding of solar flares and eruptions that 
is founded on theory and on numerical models of the behavior of magnetic 
fields from before they intrude into the solar atmosphere to their subsequent 
evolution within it. Unfortunately, our modeling abilities in these respects are 
limited by the enormous gradients and scale ranges involved (Sect. 5). For 
one thing, simulations of the emergence of magnetic flux from deep within the 
convective envelope (Sect. 10) have made substantial progress in recent years, 
but they all require twists of flux ropes that are not obviously compatible 
with the bulk of the regions observed to emerge through the photosphere: 
most solar regions emerge as bipolar regions with fibril and loop connections 
through chromosphere and corona more or less parallel to the line connecting 
the centers of gravity of the two polarities, whereas simulations would have the 
initial fields come up very nearly at right angles to the flux rope's subsurface 
direction (Sect. 10). 

Even when the field is in the corona, we have to be very careful in drawing 
conclusions from modeled field topologies and derived energy estimates. The 
main reasons for that problem are (1) the forces acting on the field within the 
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photosphere, (2) the uncertainties on vector-field measurements, particularly 
on the transverse component, and (3) the large domain that needs to be mod- 
eled to capture the connections of an active region to its surroundings (see 
Sects. Eland 8). Observationally, pre-to-post flare changes in the photosphere 
are detectable but not very strong even for the largest flares (Sect. 8.1), while 
interpreting changes in the corona suffers from the fact that it is difficult to 
disentangle coronal field changes from changed thermal conditions that cause 
different field regions to show up (Sect. 8.2). 

In view of these problems, most of our knowledge about what powers and 
triggers flares comes from observations (supported by generic models and the- 
ory; Sects. 7, 10, and 11). It appears that flares and eruptions are intrinsic to 
the emergence of magnetic field interacting with the overlying field both within 
and above the active region core field (Sects. 7 and 9). Hence, it is crucial that 
we extend our thinking from instantaneous, 2D observations and sketches to 
an intrinsically-dynamic, 3D view to capture the evolution of twisted field 
as it emerges through the photosphere, already carrying current and helic- 
ity well before doing so. Such thinking - clearly included in the sketches in 
Figs. [U HI and El and further discussed in Sect. 10 - has to elucidate why the 
large-scale polarity-inversion lines of pre-existing photospheric flux are such 
preferred sites for subsequent flux-rope emergence, which probably requires 
understanding of the deep active-region field configuration. 

Flares, flux-rope eruptions (with or without chromospheric material show- 
ing as filaments), and CMEs associated with active- region events are largely- 
overlapping populations of closely-related features of field destabilization (Sect. 3). 
It appears that the factors that determine which of these features dominates 
the evolution are matters of scale, available energy, twist, and the surrounding 
magnetic field (Sect. 9). Large statistical studies are needed to better elucidate 
the relationships between these features, and to separate the causally related 
parameters from the many others that exhibit substantial correlations. 

Essential to the occurrence of M- or X-class flares is a strong-gradient 
polarity-inversion line (SPIL) in the strong-field interior of an active region. 
I have found no cases of major flares in the absence of such a SPIL, or when 
there is not enough flux associated with it (Sect. [7]). Such SPILs are generally 
associated with flux emergence, and for the cases for which vector-magnetic 
data are available, that emergence pattern looks like a wound flux-rope that 
breaches the surface (see the well-studied example in the bottom panel of 
Fig. [Q. Such systems either flare (often repeatedly, see Sects. 4 and 6) or re- 
lax back to a near-potential state within a day or so after flux emergence ends 
(Sect. 9). 

The emergence of these flux ropes requires that most of the mass that they 
carry is drained back into the depths below the photosphere. For a spiraling, 
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undulating field this presents problems for the field-line segments that thread 
the corona on two sides with a sub-photospheric segment in between: these 
undulating field lines (Sect. [TT]) form anchors to the rising flux ropes because 
the slumping mass remains tied to the magnetic field. The pinching-off of the 
mass-carrying pockets by near-photospheric reconnection may be one of the 
mechanisms involved in the destabilization of the field configuration. 



Even as flux-rope emergence appears fundamental to the driving of ma- 
jor field destabilizations, it does not necessarily do so. Perha ps the available 
energ y can be dissipated in a series of smaller flares (e.g., ISchrijver et al 



20051 ) - although the flare magnitude spectrum per region appears statisti- 



cally to be the same (but see Sect. 4) - or by enhanced coronal heating (e.g., 



Falconer et al. 19971 ) - although the coronal energy losses appear to be inde- 



pendent of the pro perties of the vector-magnetic field in the photosphere (e.g., 



Fisher et al.lll998l ). Or, perhaps this happens because the twist of emerging 



flux ropes must be aligned with the direction of the overlying field to build up 
energy; if largely anti-parallel, reconnection may happen immediately w ithout 
the buildup of much energy (cf., simulations by lFan and Gibson! 120041 ). More 
comprehensive multi-day studies of active regions are critical to understand 
this. 



The finding that the flare-magnitude power-law spectrum of individual ac- 
tive regions is similar to the ensemble's average distribution (see Sect. 4) 
suggests that the energy-release process is the result of some process of self- 
organization of the coronal field. But it may also be a consequence of the 
sub-photospheric flux-rope shredding (Sect. 10). Repeated flaring of a region 
with continued, intermittent flux emergence supports the latter, while observa- 
tions that a flare or eruption releases only part of the region's non-potentiality 
supports the former. Understanding the consequences of these processes on the 
possible existence of a deterministic flare forecast requires systematic large- 
sample studies to differentiate the roles of the sub-photospheric and atmo- 
spheric processes. 

With respect to flare forecasting (see Sect. H21) . I note that (1) if there is 
no sign of substantial flux ropes within the photosphere for about a day, no 
major flares will occur within the next few hours, but (2) flux rope-emergence 
may lead to major flares within about a quarter of a day, and (3) may be 
associated with one or more major flares, or with a series of smaller flares, 
or with gradual dissipation of the flux-rope's energy. At present, there are 
no flare- forecast metrics with an outstanding 'skill score', although the total 
flux near a SPIL appears to be a good indicator of the flare potential of an 
active region. It remains to be seen whether the limited skill score of forecasts 
can be improved with additional observations, or whether the flare process is 
intrinsically probabilistic with the possibility of a wide range of responses to 
any given boundary perturbation. 
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2.2 A scenario for initiation of major flares and eruptions 

The literature on flares and field destabilizations, on properties of flux emer- 
gence, and on field topology and reconnection, suggests the following scenario 
for the majority of the large active-region flares that are the focus of this 
review, and whose component elements are discussed in subsequent sections. 
The top panel of Fig. 1 illustrates the components of this scenario. 

First, a large active region emerges, itself already carrying substantial net 
twist, helicity, and electrical currents, as witnessed by surface motions such 
as rotating sunspots, by its net kinetic helicity, and by the generally complex 
shape of the magnetic field. 

Subsequently, a current- carrying bundle of magnetic flux (either a single 
flux rope, or an ensemble of such ropes, perhaps formed by pre-emergence 
fragmentation) breaches the solar surface (see Fig. [T]-top). This rope may 
be a separate rise of flux through the convective envelope, or may be a late 
emergence of strands of flux from the same parent flux bundle that formed the 
active region. One flux rope, or a series of them, may emerge, forming a strong- 
gradient polarity-inversion line (SPIL) characterized by two parallel opposite- 
polarity ridges of line-of-sight flux in close proximity. At high resolution, the 
vector magnetic field across the SPIL will be highly sheared, and rotating as 
the flux bundle moves through the photosphere. 

Interaction with the near-surface convection, radiative cooling in the strati- 
fied atmosphere, and the draining of heavy plasma into surface-crossing undu- 
lations or windings of the field appear to result in a number of 'anchor' points 
for the emerging flux rope(s). Near-surface reconnection ("tether cutting") 
enables the rope to break any ties to the interior between its two extremities, 
so that part or all of it can rise relatively freely into the overlying corona. 

Models lead to the hypothesis that a rope may be unstable, and accelerate 
into an eruption, often associated with flaring, if the strength of the pre- 
existing field decreases rapidly enough with height (Sect. 9). Part of all of 
the emerged strands may participate in this, depending on the remaining sub- 
photospheric anchors. If the field gradient slows with height or changes sign, 
the flux rope may stabilize again, resulting in a confined flare or failed flux- 
rope (and, often, filament) eruption. If the field gradient is large (particularly if 
the field goes through a null, as in the "breakout" concept), the eruption may 
proceed higher, if not in fact into the heliosphere as a CME. If the winding of 
the emerging flux rope is too strong, then (confined or ejective) kinking may 
happen after any photospheric anchors between the end points are released; 
such kinking may, but needs not, transition to an ejective eruption. 

Subsequent emergence of a series of flux rope fragments is likely respon- 
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Fig. 1. Top: Schematic representation for an emerging field configuration involved 
in flaring or flux-rope eruption as discussed in Sect. 12.21 Bottom: NLFF field model 
over an active region prior to an X3.4 flare (case 8 in Table 1), rendered from a 
perspective 50 degrees from the normal view. The two main spots are connected by 
a flux rope with strong electrical currents (shown by a grey isosurface) within which 
field lines spiral from end to end on the photosphere. The Hinode vector-magnetic 
data and the extensive NLFFF modeling are discussed by Schrijver et al. (2008). 
See Sect. 5 for a discussion of NLFFF modeling capabilities. 
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sible for repeating and homologous flares, which can manifest themselves as 
leaving non-potential energy (or helicity, or electrical currents, or filaments, 
. . . ) behind in the corona if these are not (yet) unstable or released from their 
photospheric anchors. Or they may manifest themselves as repeating flares 
when subsequent strands of the flux rope emerge. 

How general is this scenario? The following sections address various aspects 
of this, but let me point out here, in summary, that the vast majority of major 
flares are associated with flux emergence (within the preceding day or so) in the 
form of tightly wound flux ropes in active regions with large-scale twist, that 
field gradients above the zone near the polarity-inversion line help differentiate 
confined from ejective filament eruptions, while the photospheric differences 
between immediately before and after the flare are generally so subtle they are 
detected only for the largest flares that have been observed at high cadence and 
then only after careful inspection (Sect. 8.1). Conversely, no major flares occur 
in regions that have a largely potential appearance, or more than about a day 
after the emergence of substantial flux ropes (Sect. 9). It remains unclear what 
role tether-cutting or anchor-cutting play in near-photospheric reconnection 
(Sect. 11). These processes have both observational and theoretical support, 
but it remains unclear how instrumental they are in the triggering of flares, 
either by themselves or in combination with the (evolving) field gradient above 
the polarity- inversion line (Sect. 9). 



3 Flares and eruptions: confined, failed, ejective 



Flares, filament eruptions, and coronal mass ejections (CMEs) are three 
aspects that manifest themselves following the destabilization of the coronal 
magnetic field. Their precise relationship remains an area of intense study, 
which is aided by the ever-increasing observational coverage of the Sun and 
inner heliosphere, yet hampered by the fact that different observatories are 
involved that generally do not have overlapping fields of view, and in many 
cases do not have continuous coverage of a given target active region. Another 
problem is, of course, that the sources of CMEs that originate on the far side of 
the Sun cannot be observed, while there is ambiguity in differentiating front- 
side CMEs from far-side CMEs - at least there was until the launch of the 
STEREO spacecraft late in 2006. 

In addition to flares and CMEs, there are the 'failed eruptions' in which 
filaments begin to rise, but th en stop their progress while twisting and writhing 



(see, e.g., iGilbert et al.l 120071 . for a discussion of this phenomenon and further 
references) within the confines of what appears to be a domain of connectivity 
for the field bounded by a separatrix; some particularly well-observed examples 
of this with TRACE are cases 1 — 7 in Table [TJ These failed eruptions are 
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Table 1 

Sample flaring regions and flare events discussed in the text. 



No. 


Time, flare magnitude 


Notes 


Sect. 


1 


2003/05/02 02:47 M1.0 


Apparently failed filament eruption obs. by TRACE/171 A. 


El 


2 


2002/05/27 18:00 M2.0 


Apparently failed filament eruption obs. by TRACE/195A. 




3 


2001/08/01 18:30 — 


Apparently failed filament eruption obs. by TRACE/17lA. 


3 


1 


2001/04/15 21:59 C5.1 


Apparently failed filament eruption obs. by TRACE/171 A. 


3 


5 


2000/11/05 00:17 C5.4 


Apparently failed filament eruption obs. by TRACE/171 A. 




6 


1999/10/20 05:53 Ml. 7 


Apparently failed filament eruption obs. by TRACE/171 A. 


3 


7 


1999/05/31 22:10 - 


Apparently failed filament eruption obs. by TRACE/17lA. 


3 


8 


2006/12/13 02:40 X3.4 


Extensive NLFFF modeling by Schrijver et al. (2008). 


El 


9 


2008/08/25 16:23, X5.3 


Ribbons distant from main P1L within AR complex. 


El 


10 


2004/07/16 13:49, X3.6 


Ribbons distant from main PIL within AR complex. 


'5: 


11 


2001/04/10 05:06, X2.3 


Ribbons distant from main PIL within AR complex. 


El 


12 


2001/10/19 00:47, XI. 6 


Ribbons distant from main PIL within AR complex. 


El 


13 


2003/03/18 11:51, XI. 5 


Ribbons distant from main PIL within AR complex. 


El 


11 


2000/06/07 15:34, XI. 2 


Ribbons distant from main PIL within AR complex. 


El 


15 


2004/07/16 10:32, XI. 1 


Ribbons distant from main PIL within AR complex. 


El 


16 


2002/08/21 05:28, XI. 


Ribbons distant from main PIL within AR complex. 


El 


17 


2004/07/17 07:51, XI. 


Ribbons distant from main PIL within AR complex. 


El 


18 


2001/06/15T10:01 M6.3 


Ribbons in adjacent AR at ~5arcmin. 


El 


19 


2002/03/14T01:38 M5.7 


Ribbons in adjacent AR at ~2arcmin. 


El 


20 


2002/07/31T01:39 Ml. 2 


Sympathetic flare in AR 4 arcmin. away. 


El 


21 


2005/09/07-17 10X, 24M 


AR 10808: 10 X-class and 24 M-class flares. 


6 


22 


2003/10/18-11/05 11X 


AR 10486: 11 X-class flares, including a record X28. 


6 


23 


2000/11/24-27 5X, 2M 


AR 9236: 5X and at least 2 M-class flares. 


m 


24 


1999/09/19 08:40 


Post-event arcade above a residual filament (TRACE/171 A) 


m 


25 


2000/07/19 23:39 


Erupting filament with decreasing twist (TRACE/ 171 A) 


1531 



See http://trace.lmsal.com/POD/bigmovies/ for a collection of TRACE movies. 

associated with flares up to low M class. At lea st some of these failed eruptions 
have no CME counterpart (e.g.. I Ji et al.ll2003l ). 



The fraction of flares that is associated with CMEs increases rapidly as 
we go from small flares to la rge X-class flares, reach ing close to 100% for 
the largest ones. For example, IWang and Zhang) (120071 ) find that 90% of 104 
studied X-class flares are associated with eruptions, and lYashiro et al.l (120051 ) 
find that the fraction of flares associated with CMEs increases from 20% for 



C3-9 flares to 100% for flares above X3. lAndrewd (120031 ) finds that 40% of 
M-class flares and all X-class flares in his sample are associated with CMEs. 

It remains unclear precisely what determines whether a large flare is asso- 
ciated with a successful ejecti o n or w ith a failed eruption. One clue is found in 
a study by IWang and Zhangl (120071 ) . who find that the 10% of X-class flares 
that they find not to be associated with CMEs tend to occur closer to the 
magnetic center of gravity of the region, leading them to suggest that erup- 
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tive flares lie closer to the periphery of the active regions. Another distinction, 
suggested by iNindos and Andrews! (12004 ). may be the field's helicity or net 
electrical current: they find that the best-fit linear force-free (LFF) fields have 
a value of the twist parameter a that is statistically larger in CME-producing 
flaring regions than in those in which the flares are not associated with CMEs, 
although there is substantial overlap in the samples. It is not clear what such 
LFF fits measure, though; a NLFFF formalism would have been appropriate, 
but that, too, is fraught with problems (Sect. 5). 



Larger (brighter and often hotter) flares tend to be associated with faster 
and /or wider CMEs (wh i ch presumably tend t o be more energet ic) ; see^Jorex^ 
ample. iKav et all (l2003h . lYashiro et ail fl2005h . lGuo et al.l(l2007h . and lGeorgouha 
(200 8|). Moreover, the fa stest CMEs originate in active regions. For example, 
Wang and ZhaneJ d2008h find that 90% of the fastest 0.5% of CMEs (with 
speeds above 1500 km/s) originate in active regions. This suggests that the 
energy in the surrounding magnetic field is involved in determining whether 
a filament eruption develops into an ejection. Other aspects are the geometry 
and gradient of the overlying field - see Sects. 0] and [9j 



4 Flare statistics from large to small 



The energy that is released during flares is converted into non-thermal 
particle populations, thermal populations, and bulk kinetic motion in case 
of ejecta. Although the initial release partitions the energy over the ther- 



mal and non-thermal partic le populations (see, e.g.. iKrucker and Lin! 12005 



Saint-Hilaire and Benzll2005l . and references therein), much of the non-thermal 
population is eventually converted into thermal energy and into the radiation 
that we see as the X-ray flare. Measurements of the radiative losses during 
flares (using, e.g., peak fluxes or fluences) thus capture much of the flare's en- 
ergy release, apart from the bulk kinetic energy in ejections and the escaping 
particle populations. 

Despite this incomplete access to the total energy involved in flares and 
eruptions, many studies suggest that peak fluxes or fluences are among the 
useful metrics for flaring activity. Distribution functions of estimated energies 
of solar flares can be approximated by power laws rather well (see, e.g., the 



overview in Table I in Aschwanden et al. 



1998, who review flare statistics from 



hard X-rays to radio as well as energeti c particle event s; see also lSu et al.ll2006l . 
for an analysis of RHESSI flares; and [Hudson! 120071 . for a compilation of 32 
years of GOES 1 — 8A flares). Figure [2] summarizes the flare-energy power 
laws from nano- to X-class flares (with energies of (9(10 31-32 ) ergs). 



The existence of power laws in not limited to flares. iRobbrechtl (120071 ) finds 
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Flare energy E [erg] 



Fig. 2. Distribution functions for estimated flare energies derived from observations 
in the extreme ultraviolet, soft X-rays, and hard X-rays. The energy range in this 
diagram spans eight orders of magnitude. From Aschwanden et al. (2000). 

essentially a power-law distribution for the widths of CMEs. Fla res with or 
without CMEs may be statistically different: lYashiro et al.l (120061 ). for exam- 
ple, suggest that power-law distributions for peak fluxes, fluences, and dura- 
tion are significantly steeper for f lares withou t CME s than for flares associated 
with CMEs. On the other hand, IWheatland reports that the waiting- 

time distribution for CMEs is like that of a time-dependent Poisson process, 
similar to that found for flares observed in the GOES records. 



In revie wing the evidence for the co mmon occurrence of power laws for flare 
properties, ICharbonneau et al.l (120011 ) conclude that their existence "suggests 
that the flaring process is intrinsic to coronal magnetic fields, even though 
the flaring rate may be controlled by extrinsic factors, such as magnetic flux 
emergence in the photosphere." The inference that the power-law distribution 
of flare properties may reflect some intrinsic property of the corona is consis- 
tent with the finding that the same power-law dis tributions a r e foun d when 
looking at subsets of flares from particular regions: IWheatlandl (l2000bl ) shows 
that the distributions of peak flare strengths for individual active regions is 
consistent with the joint pdf for all regions, and concludes that, apart from 
the frequency, 'the statistics of flares appear to be independent of the physical 
characteristics of the origi nating active r e gions' , at least within the statistical 
uncertainties. A study by iKucera et al.l (119971 ) supported an upper-limit to 
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flare magnitudes dependent on the region's spot coverage (cf., Sect. 7). They 
also found a reduced number of small flares in regions with high spot coverage, 
but sugges t that an observatio nal bias in their data requires further study on 



that issue; IWheatlandl (j2000bl ) finds no significant evidence for such an effect. 



Hudson! (120071 ) comments on the predictability of X-class flares from a cycle 
perspective. He points out that the frequency of these events through three so- 
lar cycles does not follow the Poisson statistics based on metrics for the sunspot 
cycle that appear applicable to lesser flares: there are more major flares in the 
cycle-minimum phases than there ought to be if their occurrence scaled simply 
with active-region frequency. The flaring frequency appears rather specific to 
active region properties and their sub-photospheric sources, and these appear 
to change with the cycle in ways that remain to be uncovered. 

The statistical commonality of the shape (apart from the overall frequency) 
of the flare-energy spectrum for all regions has led to the introduction of 
the concept of the ' avalanche model' or 'self- organized criticality' (SOC; see, 



e.g., the review by ICharbonneau et al.l 120011 ) in which the size of a flare is 



only statistically determined, and in which the flare's energy is always much 
less than the available total energy (thus allowing repeating and homologous 
flares with little change in the boundary field, see Sect. [5]). The SOC of the 
coronal field may develop as a consequence of the processes that inject the 
non-potential energy into the active-region coronae - discussed further in the 
following sections. We have to wonder, though, whether the 'self-organization' 
reflects a coronal property or is rather a property of the subsurface processes 
that supply the energy in the first place (see Sect. 10). With that in mind, let 
me now proceed with a discussion of what we know - or rather what we do 
not know - about flare energetics based on field measurements. 



5 Field models and topology, free energy, and helicity 



Models of solar flares would have free energy in excess of a lower-energy field 
to which it can relax, (relative) helicity, and field topology play important roles 
in determining whether a flare can occur, where it will begin its development, 
and what the final state of the coronal magnetic field can be. In order to 
derive such properties, ideally one would like to model the coronal field in 
detail based on the observable boundary conditions in the low atmosphere. 

The electrical currents in the solar atmospheres above flare-prone active 
regions are substantial, and we know that such currents can lead to significant 



chang es in field geometry (see Sect. [9]) and topology (e.g. jHudson and Wheatland 



19991 ). Unfortunately, modeling the coron al field has turned out to be exceed 



ingly difficult. Problems abound (see, e.g.. lKlimchuk et al. Il992l iMcClymont et al 
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19971 . ISchriiver et all2006al. IMetcalf et al1l2006l . iMetcalf et aDl2008l : ISchriiver et~aL 



20081 . and iDeRosa et al.ll2009l ): nonlinear force-free field (NLFFF) models or 



virial-theorem applications based on photospheric vector fields (that are sub- 
ject to problems of their own, including inversion problems and an intrinsic 
180° ambiguity for the transverse field) have to be developed further if they are 
to provide reliable es timates of the free energy, relative helicity, or topology 



Schrijver et al.l (120081 ) illustrate this very clearly. They compare field models 
for AR 10930 (case no. 8 in Tabled]) based on high-quality Hinode/SP vector- 
magnetic measurements, using a variety of NLFFF algorithms and a range 
of boundary conditions. Only one of these solutions approximates the ob- 
served corona l config uration (bottom panel in Fig. 1). A subsequent study by 



DeRosa et aJj (120091 ) based on Hinode and STEREO observations concludes 



that there is, as yet, no successful strategy for NLFFF modeling that can be 
applied under general circumstances to yield significant energy estimates and 
reliable field geometries. Consequently, in my opinion, we have yet to suc- 
ceed to a level that makes comparison of observation-based field models with 
theoretical concepts valuable to differentiate between different flare scenarios. 



When studying the changes in field properties from before to after a flare, an 
added problem to the one discussed above is that the volumes that can be ef- 
fectively modeled within computational limitations are so small that they have 
no real solar counterparts that may be considered to be closed to any physical 
quantity. Hence, the apparent convenience of having a conserved quantity - the 
relative helicity - can in reality not be exploited effectively, even if it could be 
reliably measured. 

Then there is, of course, the problem that the coronal field is anchored in 
the photosphere (which shows only subtle lasting changes during the course of 
a flare, see Sect. 18.11) . which limits the field's relaxation towards a LFF field. 
Moreover, in any given flare, only part of a region's helicity and energy may 
be involved (see Sect. 7). Consequently, even if we were able to measure the 
'free energy' in excess of either a potential field or a LFF field (at the same 
or any other relative helicity), we would only be able to place an upper limit 
on the magnitude of a flare. 



Those studies that do try to measure energy or helicity in spite of the above 
problems, typically find that the energy or relative helicity in the region suf- 
fices to power a flare or CME, and that they typically decrease from before to 
after the flare, even when measuring relative to a LFF field at the same helic- 



ity (e.g.. iBlevbel et al.ll2002l : IMetcalf et al.ll2005l: iRegnier and Canfield 



Regnier and Priestl l2007al [also iRegnier and Priestl l2007bl ] ; iBobra et al. 



2006 



2008 



Liul l2007h . Perhaps we can have more confidence in relative measurements 



than in absolute measurements? Or perhaps we are establishing that flares re- 
lease only a small fraction of the total energy associated with the field, so that 
we should not be surprised that the energy in a (N)LFF field model exceeds 
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the potential or LFF-field energy by enough to let us conclude that there is 
enough energy to power the flare. That would be a rather more spurious result 
than one would like to have. 



The application of the virial theore m to the lower-boundar y field suffers 
from its own problems. For example, I Wheatland and MetcalJ (120061 ) argue 
that the virial estimates of the free energy in active region fields based on 
the surface vector field can be somewhat improved over the simple boundary 
where that ignores Lorentz forces, but that the results are subject to very 
substantial uncertainties. An application of the virial theorem to a model field 
in the detailed study by iMetcalf et al.l (120081 ). for example, shows that the 
magnitude of the uncertainty in the virial-theorem estimate of free energy is 
comparable to that in NLFF models. 



If we cannot, yet, rely on measurements based on field models, can we at 
least rely on the overall geometry of the field, and in particular properties of its 
skeleton, such as (quasi-)separatrices? For nearly potential fields, this should 
be possible, but the largest flares and CMEs are associated with decidedly 
non-potential fields (se e Sect. ED, which unfortunately cannot be successfully 
modeled. For example, iHudson and Wheatland! (119991 ) point out that the ob- 



servation that flares occur primarily at separatrices (or quasi-separatrix layers) 
is largely unproven because of the mapping ambiguities in 3D field models. 
Where such agreements are found, we should probably even question whether 
that is because in such cases the non-potentiality (or force-free parameter a) is 
rather small, so that we are largely looking at the properties of a potential-field 
configuration. 



One property of flares that has proven quite telling as far as the reconnected 
field involved is the phenomenon of flare ribbons. They have long been asso- 
ciated with impact sites o f energetic partic les caused by reconnection (for one 
detailed study on this, see lAsai et al.ll2003l ). Generally, flare ribbons consist of 
partially- or fully-outlined fron ts moving away on either si de o f a pronounced 



polar ity inversion line; see, e.g. . Fletcher and Hudsonll200ll . and lFletcher et al. 



20051 . for discussions of particular events, with references to other studies). 
They are interpreted as the signatures of reconnection following the erup- 
tion of field through an arcade-like overlying field, leading to particle accel- 
eration and precipitation, energy thermalization, and thus eventually to the 
two- ribbon flare observable in visible light, UV, and EUV. In the case of a 
true ejection of a (frequently filament-carrying) flux rope into interplanetary 
space, one would expect that the erupting field has to break through more 
than the lowest domain of connectivity for the field, and that consequently 
more 'ribbons' should show up than only the innermost ones (which is likely 
to be the brightest and most readily observable one). It seems that this phe- 
nomenon of more distant particle impact sites has received relatively little 
attention, perhaps because of the problem of limited fields of view, perhaps 
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because their interpretation requires a large-scale field model that is generally 
not available. Nevertheless, a study of such events (of which I list some ex- 
amples in Table [TJ cases 9 — 19) may prove quite telling on large-scale field 
topology and its changes, and should be much easier with the full-Sun field of 
view of the future Solar Dynamics Observatory. Studies of such larger areas 
may also shed light on the phenomenon of the sympathetic flare (an example 
of which is case no. 20 in Table [1]). 



Li et al.l (119971 ) find that hard X-ray flare footpoints occur preferentially 
at the edges of channels of strong vertical currents as inferred from vector- 
magnetograms. This is suggestive of reconnection between a current-carrying 
rope and its surrounding field at flare onset (cf., Sect. 8.1). Higher-resolution 
Hinode data are well suited for a follow-up of that work, once the cycle picks 
up strength. 

Initial flare brightenings, whether in ribbons or in more compact impact 
sites or compact flari ng sites, genera lly lie close to the primary, strongest 



polarity-inversion line: ISchrijverl (120071 ) find that the average distance D of the 



brightest TRACE (E)UV location at flare onset (the first brightening within 
about 1 ksec from the GOES flare start time, often shown by the EUV diffrac- 
tion pattern associated with the filter support grid) to the nearest point on a 
SPIL within the region is 20 Mm. The distribution of D is strongly peaked at 
small distances, such that about 70% of the flares have D < 15 Mm, compared 
to a typical scale for the active regions of at least 150 Mm. Such distance mea- 
surements can be used i n combination with detailed models, as in the study by 



Gary and Moord (120041 ) . who discuss a filament eruption in which the initial 
brightenings lie relatively far enough from the polarity-inversion lines that 
they argue this is evidence for the initial reconnections to lie above, rather 
than below the erupting filament. I return to the role of SPILs in Sect. [7] 
and to the merits of 'breakout' and 'tether-cutting' concepts involved in the 
argument of Gary and Moore in Sect. [9j 



6 Homology and repeaters 



Among the potentially telling properties of flares we find the phenomenon 
of repeating flares for a given active region, and in particular the sub-category 
of the homologous flare. In homologous flares very nearly the same site and 
flare ribbons are involved in sequences of two or more flares. Among the well- 
studied regions of repeated flaring, we find, for example, AR 9236, AR 10808 
(among the most flare-productive regions recorded based on the GOES soft 
X-ray data), and AR 10486 (see entries 21 — 23 in Tabled]). 

The repeating (possibly homologous) flare suggests that either not all of 
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the available energy is released in a flare, or that energy continues to be sup- 
plied to the region so that a series of flares of various magnitu des can occur. 



Evide nce is found for both of these possibilities. For example, iBleybel et al. 



( 120021 ) study a case in which the post-flare configuration remains inconsistent 
with a LFFF model (either based on the pre-event or the post-event rela- 
tive helicity as determined by a Grad- Rubin based NLFFF code). Multiple 
studies have reported on contin ued flux emergence between successive ho- 



mologous flares (see, for e xample. iRanns et al.l l2000l: iNitta and Hudson! 12001 



Sterling and Moorell200ll ; iTakasaki et al.l 12004 ; and lDun et al.l 120071 ). Some of 
these repeater events clearly show that not all the available energy is taken 
away in an event. One such example is that of case 24 in Table [1] which shows 
a post-eruption arcade over a residual filament structure. 

Some studies argue that systematic motions, such as the persistent rotation 
of sunspots, are invo lved in rebui l ding t he energy for homologous flares (see 
als o Sect.fTUl). such aslZhang et al.l (120081 ) (see also the related numerical model 
by DeVore and Antiochoj 20081 ). The phenomenon of the rotating sunspot is 
generally associated with ongoing flux emergence, how ever (such as in the case 
of the X17.2 flare Oct. 28, 2003, one of the subjects of Izhang et al1l2008l . part 
of entry 22 in Table 1), and without flux emergence, major flares rarely, if at 
all, occur (see Sect. 8). Besides, it is hard to estimate whether flux em ergence 



or spot rotation adds most energy (e.g., iRegnier and Canfieldl 120061 ) in view 
of the discussion in Sect . [5j The reformatio n of loop systems involve d in long- 
distance, even trans-equatorial connections (IKhan and Hudsonll2000l ) is harder 
to interpret, but perhaps also hinting that not all energy is removed during 
eruptive events. 



7 Patterns in the surface field 



The most frequently used observable of the magnetic field involved in im- 
pulsive and eruptive events is the surface line-of-sight field, largely because it 
is readily measured (and has been, routinely, for a long time). Increasingly, 
vector-magnetic measurements have become available. Often, the frequency of 
either of these measurements is rather low, in part because of observational 
limitations for ground-based observatories. Consequently, many studies have 
focused on finding properties of patterns of the surface magnetic field in snap- 
shot line-of-sight magnetograms that might shed light on the powering and 
triggering of impulsive phenomena. 

One type of study focuses on statistical properties of the magnetic field of 
active regions. First, it appears that M- and X-class flares occur in sunspot- 



containing regions with ver y few e xceptions. The work of, e.g.. lDodson and Hedeman 



(Il970l ) and iRuzdjak et al.l (119891 ) points out that whereas flares do occur in 
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Fig. 3. Scatter diagram of peak flare flux densities (GOES 1 — 8 A passband in 
W/m 2 ) versus the absolute total flux near the strong-field, high-gradient polarity 
inversion, R (in relative units). The histograms in this figure show the 25-th, 50-th, 
and 75-th percentiles of the flare distribution as function of R. The dotted straight 
line at -F pe ak = 1-2 x 10 -8 R has 99% of the observed flares below it. The only flare 
that clearly lies above this line is an X6.2 flare that started around 2001/12/13 
14:30 UT. The value of log (R) in this active region decreased from 4.2 half a day 
prior to the flare to 3.7 at the time of the flare, but at 4.2 would still lie well above 
the dotted line. This region has a narrow fibril of one polarity embedded within 
an area of strong field of the opposite polarity, which likely causes the measure R 
to strongly underestimate the true flux involved (see footnote). This may (in part) 
explain its unusual position in this diagram. From Schrijver (2007). 

spotless regions, these account for only about 2% of all flares while the largest 
documented flare in soft X-rays of this type was classified as Ml. 2. This is the 
only documentated case I found for which a major flare in soft X-rays occurred 
in a spotles s region; a l l of th e 289 M- and X-class flares, for example, that were 
studied by ISchrijverl (120071 ) - which were the basis of the study resulting in 
Fig. 3 - have spots associated with them. 



Another statistical property was pointed by by lCanfield and Russelll ( 120071 ) . 
who find that a tessellation of active-region magnetograms at a few arc-second 
resolution yields log-normal distributions, suggesting that the flux bundles 
were subjected to repeated random fragm entation during their rise to the sur- 
face (compare the study by IBogdanlll985l . on a similar size distribution of spot 
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areas). They find no significant correlation between the degree of fragmenta- 
tion and either the flare rate of the active regions or the mean density of kinetic 
helicity, concluding that it is more likely that large-scale helicity in the flow is 
related to flaring than sm all-scale inte r action s between convection and rising 



flux. On the other hand, lAbramenkol (120051 ) finds that the indices of power 



spectra of the photospheric magnetic field of flare-prone regions increase from 
— 1.7 for flare-quiet up to as steep as —2.3 for regions producing high X-class 
flares (see her Fig. 5). It is not clear how the preceding two results are to be 
reconciled. 

Some of the salient features in the active-region magnetic field have received 
far more attention than the above-mentioned statistical properties, in partic- 
ular the properties of the polarity-inversion lines. Despite all this attention, 
no measure (either single or combined) has been identified as being partic- 
ularly well correl ated with the flare activity of active regions. For example, 
Georg ouiisl fl2008h studies 23 eruptive active regions and finds that for active 



regions with 'intense polarity inversion lines' (as measured by a metric that 
combines fluxes and field-line lengths from connectivity matrices between pos- 
itive and negative fluxes in magnetograms ) tend to have both stronger flares 



and faster CMEs. IWang and Zhang! (120081 ) find that the length and number 



of the polarity-inversio n lines (PILs) a re go od indicators of fast vs. slower 



CMEs. The findings by iFalconer et al.l (120061 ) (expanding on Falconer et al. 



2002, 2003) suggest a more complex dependency: for a sample of active re- 
gions with one well-defined, dominant polarity-inversion line they find that 
several measures are correlated similarly to the region's flare potential. These 
measures include the length of the strong-shear segment of the PIL, the length 
of the main PIL with a strong field gradient across it, the total net current, 
and the best-fit a from a LFF field based on a comparison of the photospheric 
vector field with the model field footprint. They suggest that the best metric 
for the total magnetic free energy is, in fact, a combination of total flux and 
total twist. 

The ambiguity in findings in general is likely in part the result of the cor- 
relation between multiple parame ters, and disentangling these requires large 



statistical samples. The study by iLeka and Barnes! (120071 ) is an excellent ex 



ample of that, but still leaves us with an unclear picture. They analyze 496 
active regions to establish which parameters of the photospheric field are most 
successful in differentiating flaring from non-flaring regions. They measure 28 
different properties from the vector-magnetic field measurements and differ- 
entiate flaring from non-flaring regions within 24 h of these magnetograph 
observations. Their discriminant analysis has a success rate of 93% for Ml 
flares or larger (compared to 91% when simply labeling all regions as flare- 
quiet, see Sect. fT2l) . or 80% for CI or larger (compared to 70% in the absence 
of any information on the region). They find that the best-performing single 
parameter for large flares is a proxy for the magnetic free energy in the region; 
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Table 2 

Likelihood of X or M flares within 24 h of the determination of the unsigned mag- 
netic flux 1Z within about 15 Mm of high-gradient, strong-field polarity-separation 
lines (see Schrijver, 2007, for a detailed description and normalization factor). The 
likelihoods are sampled at 0.5 dex intervals, but vary smoothly with 1Z. Also listed 
is the maximum expected flare class. 
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other top performers are variables that measure the field's various properties 
integrated over the entire active region. 



Another example of a statistical study is that by lSchrijverl (120071 ). who looks 
at the magnetic properties of regions associated with almost 300 M- and X- 
class flares, and compares that with the properties of 2,500 randomly selected 
active regions, all within 45° of disk center. He notes that the regions with 
large flares all have a pronounced polarity-inversion line with a strong gradient 
across it. He subsequently measures the unsigned flux 1Z within about 15 Mm 
of such strong PILs (SPILs), and this is a good indicator for flares that allows 
prediction of the largest flare to be expected from the region (see Fig. 3 and 
Table 2; I return to the issue of predictability and the skill of various metrics 
in Sect. H21). 



Based on the inspection of many full-disk S OHO -MDI magnetogram se- 
quences for regions with large flares, ISchrijverl (120071 ) argues that the SPILs 
appear to be associated with emerging flux, in particul ar with emerging flux 



ropes (also supported by the det ailed field modeling by lSchrijver et al.l 12008 



cf. Fig. 1). IWelsch and Lil (120081 ) show that there is a pronounced tendency 
of increasing 1Z to be correlated with increasing total unsigned flux, i.e., with 
signs of flux emergence. Given the resolution of full-disk MDI magnetograms 
and the likelihood that emerging flux immediately cancels partly against pre- 
existing flux, this is strong evidence for flux emergence as driver for the SPIL 
flux metric 1Z. More evidence for the association with flux ropes in particular 
is found when looking at shear at such SPILs. 



Shear at the polarity-inversion line is ofte n associated wit h flaring activ- 
ity (examples are a statistical study by , e.g.. ICui et al.1 120071 and individual 



event studies b y Kurokawa et al. 2002 . Brooks et al. 20031 . Dun et al. 2007 



Schrijver et al. 2008 . Okamoto et al. 2008 . and Lites 20081 ; also coronal field 



studie s described in Sect. |9]). The same holds for eruptions: IVrsnak et al. 



( 119911 ) report that the pitch angle of Ha structures outlining filaments is telling 
regarding the eruption potential: prominences with twist angles exceeding 50° 
and heights exceeding 80% of the footpoint separation erupted, while those 
with angles below 35° and heights below 60% of the footpoint separation were 
stable. 



Manchester! (120071 ). based on simulations of an emerging flux rope, argues 
that shear motion at the polarity-inversion line in emerging flux is a natural 
consequence of the Lorentz forces that develop as a rope enters the corona. 
This may explain the general correspondence of field shear and shear flows. 



8 Flare-related changes in the magnetic field 



8. 1 Before and after in the photosphere 



The conversion of large-scale electromagnetic energy in the process of a 
flare or filament eruption changes the field configuration over active regions, 
and this should in turn have consequences for the observable vector-magnetic 
field in the photosphere. Until recently, these changes proved elusive (other 
than temporary changes in line profiles during the flares, which are mostly if 
not entirely a consequence of the energy deposition in the near-atmospheric 
laye rs by populations of e nergetic particles generated early in the flare, see, 



e.g., ISchrijver et al.ll2006bl ). 



Part of the problem in finding flare-related field changes is that it cannot re- 
ally be done with vector-magnetic measurements based on spectrograph scans, 
simply because these scans take too long to complete relative to the evolution 
of the photospheric field between two successive magnetograms with spacings 
that are much longer than the impulsive phase of the flare or eruption. For 



example, iDun et al.l (120071 ). using a series of 12 vector-magnetograms in the 
course of five days of AR 10486 (case 22 in Tabled]), find that magnetic shear 
angles, electric current density, and current helicity increase at the sites of 
impulsive flares at least a day prior to the main (X17 and X10) flares. But 
they attribute these changes to the emergence of twisted flux ropes rather 
than to changes associated with the flares per se. Although they say their 
metrics decrease after the main flares, the observational cadence is too low to 
unambiguously differentiate flare-related relaxation of the field or the passage 
of flux ropes through the photosphere prior to the post-flare magnetograms. 

The analysis of line-of-sight magnetograms taken at much higher frequency 
has been successful in identifying differences in the magnetic field just before 
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and after major flares. ISudol and Harveyl (120051 ). for example, find: a) that 
the time scale for changes in the l.o.s. field is less than ~ lOmin., some even 
less than one minutes, but in all cases abrupt; b) in many cases these changes 
are close to the noise but in some cases very significant; c) the changes persist 
for at least the two hours past the peak flare time studied by them, and in 
one case for at least 5 hours. They note that most of these observed field 
changes occur in spot penumbrae, and that for the three flares for which 
TRACE coverage is available, these changes occur at sites of flare ribbons in 
the (E)UV. The authors argue that these field changes are lasting, and thus 
likely not the result of atmospheric changes related to the particle impact 
or heat conduction leading to the flare-ribbon formation; perhaps the energy 
deposition occurs (very nearly) along the current- carrying field lines involved 
in a reconnection event. The fact that these field changes are often close to 
the detection threshold even for major flares will make it difficult for field 
modelers to reliably deduce changes in the field energy from before to after 
flares in general - see Sect. 5). 



Wang et al.l (119941 ) show cases in which the shear angle across the polarity- 



inversion line of active regions increases during, and persists after a set of 
five X-class flares, with a time resolution that ranges from about 4-10 min. 
in the best cases to about 45 min. They propose as one possible interpreta- 
tion that new flux emerging just before the flare could be the cause of this 
increase, but comment that it is hard to understand how that could change 
the shear over more than 10,000 km in such a short time. On the other hand, 



Chen et al.l (119941 ) report on shear studies of 18 regions around M-class flares; 



they find no detectable changes, suggesting that perhaps only the largest flares 
are associated with observable changes. 



Even very large flares and eruptions are associated with relatively weak or 
very localized changes in the surface magnetic field. This is, of course, part 
of the problem of measuring significant changes in the field energy, helicity, 
and topology around the time of such field destabilizations: too often, the field 
changes that are reported are mixed in with intrinsic changes in the photo- 
spheric field caused by flows (horizontal and, often, vertical associated with 
flux emergence) owing to the substantial time intervals between successive 
magnetograms, thus masking the changes associated with the flare or erup- 
tion only. The potential of, e.g., the 1-min. cadence magnetogram sequences 
for SOHO/MDI has not yet been utilized. Moreover, we shall need sensitive 
(vector-) magnetographs, preferably for the chromospheric field, with high ca- 
dences and sizable fields of view in order to learn about vector field changes 
at the lower boundary at times of flares and eruptions. 
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8.2 Before and after in the corona 



There do not seem to be many studies that compare the coronal appearance 
before and well after a major flare (well observed in, e.g., case 24 in Table 1). 
There are some studies, however, that look at loops early in a flar e /eruption 



and in the late stages of the post-flare arcade. Among these is lAsai et al 



(120031 ). who find that during the course of the initial stages of the flare, the 



conjugate footpoints at first indicate st rongly sheared fi eld lines, and later 



reveal less sheared coronal connections. ISu et al.l (l2007al ) study TRACE ob- 
servations of 50 X- and M-class two-ribbon flares, with well-defined, largely 
parallel ribbons. They find that 86% of these flares show a general decrease in 
the shear angle between the main polarity inversion line and pairs of conjugate 
bright ribbon kernels. They interpret this as a relaxation of the field towards 
a more potential state because of the eruption that carries helicity /current 
with it, but one can readily argue that a similar decrease in shear angle would 
be seen if sequentially higher, less-sheared post-flare loops lit up with time as 
the loops cool after reconnection. These results are consequently ambiguous: 
they may show a decrease in shear, or they may reflect that flares generally do 
not release all available energy and part of the flux-rope configuration remains 
(cf., Sect. |BJ. 



Su et al.l (l2007bl ). from a study of a sample of 31 two- ribbon flares with 
CMEs, find that the primary parameters describing the peak flare brightness 
and CME terminal velocity are the total flux in the region, and the change 
in shear angle between the early and late phases of the flare, and - to a 
significantly lesser degree - the initial shear angle. The change in shear angle 
may reflect a real decrease in shear and associated energy, but - in view of the 
preceding paragraph - may also be a measure of the residual shear after the 
eruption; hence, it remains unclear what the measured shear angle is telling 
us in this case. 

Compounding the efforts to understand the field evolution are the relative 
roles of twist and w rithe, and their coupling. One study on this topic by 



Romano et al.l (120031 ) reports that the number of turns around an erupting 
filament initially is about 5, decreasing to 1 later into the eruption (entry 25 
in Tabled]). Here we may see an exchange of twist and writhe. 



9 The overlying coronal field 



Numerical MHD models continue to shed light on the potential causes of the 
instability that leads to flares and, in particular, filament eruptions, as always 
with the clear caveat that such modeling is an approximating abstraction of 
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reality. Some of these models rely on the expected general slowness of reconnec- 
tion in a solar atmospheric plasma, except for particular field configurations, 
specifically null points. The so-called breakout model, for example, "exploits 
a vulnerability of multi-polar configurations, which consist of two or more dis- 
tinct flux systems separated by null points in the corona, to rear rangements of 
the magnetic field's connectivity" ODeVore and Antiochosll2008l ). Specifically, 
rapid reconnection is anticipated as flux is pushed into a null region. O ne re- 
cent example of such fields is discussed by lDeVore and Antiochod (120081 ). They 
simulate a multi-polar configuration with a null in the corona over a model 
active region subjected to rotational flows (somewhat like rotating sunspots). 
Their simulation run extends to cover three homologous eruptions in their 
breakout concept. These, by the way, are confined, in the sense that they 
erupt only from the active region, but not into interplanetary space. 



Another concept is the 'tether-cutting' model (originally described by lMoore and Roumeliotis 



19921 ). in which (near-)photospheric reconnection associated with flux emer- 
gence removes the 'anchors' to flux rope that can subsequently reconnect with 
the overlying field, or even erupt. This is discussed further in Sect. [TTJ 



A different possibility to destabilize a field configuration is related to the 
kink and torus instabilities. The kink instability occurs when the helicity, 
current, or w inding-number for the fi eld in a flux rope exceeds a critical value. 
For example. iFan and Gibson! (120041 ) perform line-tied MHD simulations of an 
emerging flux rope into a pre-existing configuration. For a case in which the 
twisted field in the top part of the rising rope has the same general direction as 
the overlying arcade, they find that the configuration becomes unstable after 
the coronal field reaches 1.76 full turns about the rope's axis. The erupting 
rope exchanges twist for writhe as it evolves as is often - but certainly not 
always - seen in filament eruptions. In a case in which the rope's top-most 
winding relative to the overlying arcade is reversed, no flux rope builds up 
in the corona as reconnection proceeds from the very first emergence, and 
consequently no instability develops. 



The torus instability occurs whenever the gradient in the surrounding field 
is strong enough, so that the forces exerted on a flux rope s ubject to a pertur- 
bation cannot be contained by the overl ying magnetic field. iKliem and Torok 
(120061 ) (also, e.g-. lFan and Gibsonll2007l ) discuss the torus instability, and find 
a critical field gradient beyond which the configuration is unstable. They argue 
that line tying of the field to the photospheric boundary in this case helps to 
stabilize a configuration, arguing that it may even stabilize an emerging flux 
rope entirely until it is semi-circular and thereafter may help to destabilize it. 



Torok and Klieml (120071 ) discuss a unifying approach to slow and fast CMEs 



in which the torus instability drives the CME, and in which the gradient in 
the surrounding field determines the acceleration profile and eventual velocity 
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of the ejecta: fast CMEs for rapid decrease (as typical of ARs) arid sl ow for 
gradual decrease (as often over quiet Sun). Ilsenberg and Forbes! (12007) star t 
from a similar field configuration as developed by lTitov and Demouiinl (119991 ). 
but perform an analytical analysis for eruptions without twis t (i.e., without 
the subsurface current introduced by Titov and Demoulin 19991 ). They confirm 
the instability of the configuration when the flux rope length exceeds a few 
times the diameter of the flux rope. 



Following the proposition of the torus instability, IWang and Zhang! (120071 ). 
based on potential-field models, find that 'the ratio of the magnetic flux in the 
low corona to that in the high coro na i s syste matically larger for the eruptive 



events than for the confined ones'. iLiul (120081 ). for a sample of ten flares, find 



a low field gradient for failed eruptions or kinking, and a steeper gradient for 
successful eruptions. The power-law indices for the field decrease with height 
for the failed eruptions and for the kink and torus eruptions combined are 
— 1.62±0.05 and — 1.91±0.15, respectively, for height ranges from log (/i [km]) = 
3.7 — 4.7. The use of potential- field models here is an oversimplification and 
the findings may say more about the distribution of the surrounding field on 
the photosphere than about the actual field gradient in the low corona. 



The above-mentioned studies deal with CMEs in which flux ropes that may 
contain filaments are part of the story. It remains to be seen how this can apply 
to flares, or to eruptions with failed ejections, but there is no obvious reason 
why the above arguments would not apply to compact flares as much as they 
do to more extended flux-rope cum filament systems. 



In view of the theoretical and numerical considerations described in the 
preceding subsection, it is no surprise that multiple observational studies look 
at the coronal field configuration. Of course, this brings us back to the problem 
that the coronal magnetic field configuration cannot be known really well 
with present tools (Sect. E]). This is particularly bothersome because it is 
clear that the coronal fields in flare-prone regions differ markedly from the 
potential field bas ed on the same vertica l (or line of sight) photospheric field; 
see, for example, ISchrijver et al. (l2005h . This does not mean that currents 
pervade much of the active-region corona, of course, but at least points out 
that somewhere there must be currents strong enough to warp the field relative 
to a pure potential field everywhere from the core field to the peripheral field. 
It appears that such currents survive for t ypically only a day^ while powering 
one or more flares during t heir decay (e.g.. 
20041 : ISchriiver et al.lliooih . 



Pevtsov et al.lll994t iBurnette et al 



How common is a strongly non-potential field? IBurnette et al.l (120041 ) com- 
pare photospheric vector-magnetograms and coronal configurations as seen by 
YOHKOH's SXT. They compute LFF models that best match the coronal 
geometry, and find that over 90% of the regions studied are fitted well by such 
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a model. iNindos and Andrews! (120041) find that t his is true for only 60% of 
their sample. Furthermore, iBurnette et al.l (120041 ) find that the a parameter 
derived from the vector-magnetogram correlates well with that estimated for 
the coronal configuration. Apparently, the currents that are responsible for the 
large-scale non-potentiality of the coronal field close through the photosphere, 
while small-scale structure of the a parameter in the photosphere often does 
not contribute significantly to the large-scale field configuration. 



The electrical currents involved in non- potential regions app ear to emerge 
with the flux, according to several studies. iPevtsov et al.l (120031 ). for example, 
use line-of-sight MDI magnetograms and LFF field fits to EIT EUV images 
to estimate the helicity evolution in emerging active regions. They find that 
the helicity parameter a starts essentially at zero on first emergence for all six 
regions that were well observed, and subsequently converges to a maximum 
plateau value over the following 1.5 d. They argue that this means that the 
twist is in fact injected into the corona by 'the spinning of the active region 
polarities driven by magnetic torque from below'. Note that this means on 
first eme rgence, the coronal fie l d is ap parently nearly potential, with a ~ 0. 
Work by lLongcope and Welschl (120001 ) supports this finding. 



The studies in the literature, both observational and theoretical, suggest 
that whether a flux rope erupts not at all, evolves into a "failed eruption", 
or makes it into the heliosphere as a CME depends critically on the structure 
of the overlying field, in particular on the gradient of the field with height, 
but likely also the direction of the overlying field relative to the twist of the 
emerging flux rope. 



10 Subsurface processes, emerging flux, and rotating spots 



Currents through the photosphere are often not neutralized when integrat 



ing over either of the two polar ities. It has been argued (e.g., iMelrosd 11991 
Melroseill995l . IWheatlandll2000al ) that consequently these currents cannot have 
been induced by non-current carrying flux emergence, because in that case the 
currents should balance within each polarity. Hence, instead, at least some of 
these currents must have emerged from below the photosphere. This notion 
is compatible with the observationa l evidence that the Sun's mag netic field 



is often twisted on emergence (e.g., 


Kurokawa 


1987. 


Tanaka 


1991 


[see, e.g., 


his Fig. 51, 


Leka et al. 


1996. 


van Driel-Gesztelvi et al. 


1997. 


Ishii et al. 1998. 


Lopez Fuentes et al. 


2000. 


Okamoto et al. 


2008). 


Brooks et al. 


2003) study in 



detail, with high spatio-temporal resolution, the conditions of two emerging 
flux regions, one with a small flare, and one without, in the field of view of 
the SVST. They find that the flaring region is associated with the emergence 
of what looks like a twisted, i.e., current-carrying flux rope (compare Fig. 1), 
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Fig. 4. Select examples of proposed 3D configurations to explain the motions of pores 
around a spot (top left), of the complex motions of multiple concentrations of flux 
in an emerging region (top right), and the photospheric crossing of an emerging flux 
rope (bottom two panels). In these sketches the possibility of efficient supra-photo- 
spheric reconnection (as in the Ellerman bombs of Fig. 5) has been ignored. (Top 
left: from Ishii et al., 1998; top right: Kurokawa et al., 2002; bottom: from Low, 
2001) 

while the other has no such appearance. 



Not only emerging currents contribute to a region's free energy, however. 
Also syst ematic flow patte r ns at, or just below, the surface can do so. For 
example, iNightingale et al.l (120071 ). report that rotating sunspots are associ- 
ated with almost all of the X-flares observed by TRACE since April 1998, and 
many of the M-flares. 

It is difficult to differentiate twist in the emerging flux from twisting flows. 
Support tha t these twists a re ass ociated with flows, at least below the surface, 
is found by iKomm et al.l (120051 ). who estimate the kinetic helicity density 
(v • V x v) below active regions using helioseismic measurements. They find 
that the maximum helicity density in the outer 2% of the Sun correlates 
'remarkably well' with the X-ray flaring activity, perhaps - as they point out - 
because of the link between helicity and twist in the subsurface magnetic field. 



Whereas sunspot rotation has been argued to add energy to the coronal 
field, either because of the work done by the (sub-)photospheric flows or be- 
cause it reflects the emer gence of a twisted, current-carrying flux bundle, the 
study by lSchrijverl (120071 ) shows that all X-M flaring regions contain SPILs. As 
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all X-flaring regions appear to contain rotating sunspots, it may well be that 
these are related phenomena, perhaps because many, if not all, flux bundles in 
a flare-prone region have significant electrical currents associated with them 
induced by preferential flow systems. 



There appears to be a net relative helicity in each hemisphere, negative in 
the north, positive in the south, corresponding to left-handed screws in the 
north and right-handed ones in the south: the statistical study by lPevtsov et al. 
(119951 ) shows a 2:1 preferences of the sign of he licity for the hemispheres, ay - 



eraged over all latitudes and for two cycles (see lLongcope and Pevtsovll2003l ) 



This twist, if existent already below the surface , is likely to be affected by tur 



bulent convective buffeting of rising flu x (see lLopez Fuentes et al.l |2000| . for 



a discussion of an example of this). As lLongcope and Pevtsovi (120031 ) argue, 
this would trade twist with writhe of the emerging flux tubes, without affect- 
ing relative helicity. The writhe should affect the tilt angle of dipole axes on 
emergence, if not cause kinking of the subsurface flux arches into what may 
be seen as complex, or even (5-class regions. This suggestion has been made, 
but I have not found this empirically investigated. 

On the theoretical side, a persistent problem is that even state-of-the-art 
numerical experiments still have difficulty in successfully bringing flux bun- 
dles from some depth up into the photosphere without them explodin g in 
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the s trongly stratified uppermost layers of the convective envelope (e.g 
20041 ). Note that among the stabilizing forces in this proce ss we find increasing 
curvature, the Cor i olis forces, and Reynolds num ber (e.g., Abbett et al.l 2000 . 
Abbett et ai1l200ll . iMagaralliooi ICheuneJl2006h . Models suggest that twist is 
needed for the flux bundles - thus turning them into flux ropes - to survive the 
rise. On the other hand, it appears that modelers need to impose more twist 
than warranted by observations, as most active regions emerge with archfila- 
ment systems mostly aligned with the eventual dipole axis of t he region rather 
than nearly perpendicular to it as in many simu lations (e.g. 
20051 : iMagaral [20061 : lManchesterll2007l : lFanll2008h . 



Archontis et al 



Starting with t oo m uch twist subjects the emerging flux ropes to kink 
instabilities. iFanl (120081 ) discusses models of rising flux tubes through a ro- 
tating medium stratified like a convective envelope. She shows that if these 
non-axisymmetric flux tubes rise with a twist of the observed hemispheric 
preference that is strong enough to prevent substantial shredding during the 
rise, the twist converts into a writhe that has the opposite sign than the ob- 
served preferred tilt angle. If the twist is weakened by only a factor of two 
this effect becomes weaker than that of the Coriolis force, and the net result 
is a tilt that is comparable to the observed one. But the consequence of this 
is that the rising 3D tube loses a lot of its flux by shredding as it moves buoy- 
antly through the surrounding medium. In view of this, one sees this scenario 
develop: weakly twisted tubes explode and cause network field; moderately 
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twisted regions have a trail of shredded field behind them, perhaps causing 
active-region nesting; and strongly twisted field should come up coherently, 
perhaps subjected to a kinking-instability leading to 5-spots. 



11 The role of near-photospheric reconnection 



With mounting evidence for the emergence of flux ropes as a characteristic 
ingredient in major solar flaring, it is instructive to look at some of the highest- 
resolution studies made of the process of flux emergence. These observations 
suggest that emerging magnetic flux often does not emerge as a simple arch, 
as traditionally studied, but rather as an undulating field that transits the 
photospheric layer multiple times in sea-serpent fashion. In such a geometry, 
the magnetic field reaching into the chromosphere might drain its plasma 
load into sub-photospheric dips in the field, with trapped material that would 
not allow those field segments to rise. But perhaps near-surface reconnection 
allows the heavily mass-loaded pockets of magnetic field to be pinched off 
(thus creating a variant of the tether-cutting concept in which the emerging 
and pre-existing flux ropes are, in fact, one and the same). 



Lited (120081 ) discusses such observations of emerging flux made with Hin- 



ode's SP, NFI, and BFI. He argues that the observations of opposite-polarity 
patches moving towards each other, with downflows seen in the spectral data, 
present an opportunity for plasma to drain back into the photosphere, and 
then be pinched off from the supra-photospheric configuration, which can 
then rise into the co rona. Recent numerical simulations of emerging flux (e.g., 
Cheung et al.ll2008l ) support such an interpretation. 



The concept of pinching off the heavy, sub-photospheric dips in emerg- 
ing, twisted magnetic field (such as sk etched in the lowe r panels of Fig. Hj) 



has been suggested before (also, e.g., iPariat et al.l [200J). One example of 



this is the interpretat i on of Ellerman bombs as shown in Fig. [5j taken from 
Schmieder and Pariat fl2007f ). 



12 Instability, catastrophe, and predictability 



Let me, finally, address two issues before ending this review. First a com- 
ment on the nature of the instability. We have to ask whether this instability 
can be the result of a bifurcation that occurs in certain conditions under which 
there is more than one state for the field to be in given the observed photo- 
spheric vector field. The existence of more than one field configuration for the 
same lower boundary could be the explanation for why lasting photospheric 
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Fig. 5. Conceptual illustration (by Schmieder and Pariat, 2007) of emerging flux, 
with undulating strands of flux with multiple photospheric crossings. In this concept, 
the approach of pairs of corresponding opposite-polarity patches on either side of a 
dipped field line could lead to reconnection (in what might appear as an Ellerman 
bomb, or EB) that pinches off the heavy sub-photospheric segment, thus allowing 
the flux to rise. For field configurations that can destabilize, twist and or kinking 
appears needed, as illustrated in Figs. Q] and HI 



changes from pre- to post-flare states are hard to find (although they have 
been reported on, see S ect. 8.1). Not much seems to have been written on this 
topic, but interestingly, llsenberg and Forbes! (120071) point o ut that they 'have 
found evidence that each T&D Titov and Demoulinlll999l | equilibrium has a 
corresponding toroidal equilibrium with exactly the same boundary conditions 
at the solar surface, which is unstable when the T&D equilibrium is stable and 
vice versa,' thus pointing to a catastrophe scenario. Among the other studies 
that address the existence of multip le solutions (some stable, some unstab le) I 
mention the catastrophe m odels by Ivan Ballegooijen and Martensl (Il989h and 



Forbes and Isenberg) (119911 ) . Whether the possibility that the nearly force- free 



state of the coronal field allows multiple solutions to (very nearly) the same 
boundary condition will need additional work in order to be able to assess 
these effects in the context of field instabilities. 



The second point relates to the predictability based on proposed metrics 
for fla re likelihood. A landmark study in this is the work by iBarnes and Leka 
(120081 ) who evaluate the skill scores of forecasts base d on four metrics: total 
flux, the "total excess energy" (ILeka and Barnesll2003l ). the S PIL flux metric 7Z 



(lSchriiverll2007l) . and the "effective connected magnetic field" (IGeorgoulis and Rust 



2007h . They limit their sample to regions within 30° from disk center to create 



a relatively uniform set of observing conditions. They show that the highest 
success rate for any of these single parameter, 1Z, is 92.2%, whereas a uniform 
forecast that nothing will every produce an M-class flare or larger is 90.8%. 
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They conclude that none of these parameters are robust as forecasting tools. 



13 In conclusion — looking forward 

This, then, brings us back to the issue of "understanding." If we could 
deterministically forecast flares and eruptions, then we certainly could claim 
success. But we might also have reached a deep understanding if this were 
not the case: if flux-rope emergence is the driver, or if self-organization occurs 
within the corona, then no such thing as a "deterministic forecast" based on 
present-day observables is feasible in principle. How do we find out if either or 
both of these possibilities apply? It seems that statistical studies for observa- 
tions from interior to high corona are essential. The relatively slow process of 
emergence and the rapidity of flare onset also point to multi-day, high-cadence 
coverage for such coordinated observations. 

These observations need to be complemented by numerical experiments of 
flux emergence with a realistic treatment of the interior, of the near-surface 
layers, and of the corona. Moverover, we shall need to learn how to model the 
atmospheric field in order to measure properties such as energy, helicity, and 
topology. 

Coordinated multi-instrument observations of the Sun in the era of the Solar 
Dynamics Observatory, supported by extensive spectroscopy and numerical 
experiments, hold the promise of a breakthrough for our understanding of 
magnetic instabilities involved in flares and CMEs. 

Acknowledgments. I thank Hugh Hudson, Nariaki Nitta, and Alan Title 
for very helpful comments on early drafts of this manuscript, and the referees 
for very constructive suggestions to improve the presentation and to some key 
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